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Diet is associated with cancer prognosis, including head and neck cancer (HNC), and has been hypothesized to influence 
epigenetic state by determining the availability of functional groups involved in the modification of DNA and histone 
proteins.Thegoalofthisstudywastodescribetheassociation between pretreatmentdietandHNCtumor DNA methylation. 
Information on usual pretreatment food and nutrient intake was estimated via food frequency questionnaire (FFQ) on 
49 HNC cases. Tumor DNA methylation patterns were assessed using the lllumina Goldengate Methylation Cancer Panel. 
First, a methylation score, the sum of individual hypermethylated tumor suppressor associated CpG sites, was calculated 
and associated with dietary intake of micronutrients involved in one-carbon metabolism and antioxidant activity, and 
food groups abundant in these nutrients. Second, gene specific analyses using linear modeling with empirical Bayesian 
variance estimation were conducted to identify if methylation at individual CpG sites was associated with diet. All models 
were controlled for age, sex, smoking, alcohol and HPV status. Individuals reporting in the highest quartile of folate, 
vitamin B12 and vitamin A intake, compared with those in the lowest quartile, showed significantly less tumor suppressor 
gene methylation, as did patients reporting the highest cruciferous vegetable intake. Gene specific analyses identified 
differential associations between DNA methylation and vitamin B12 and vitamin A intake when stratifying by HPV status. 
These preliminary results suggest that intake of folate, vitamin A and vitamin B12 may be associated with the tumor DNA 
methylation profile in HNC and enhance tumor suppression. 



Introduction 

Head and neck cancer (HNC) is the sixth most common type 
of cancer and accounts for approximately 5% of all new cancer 
diagnoses worldwide each year.^'^ These malignancies develop 
from the epithelial tissue of the oral cavity, nasopharynx, oro- 
pharynx or larynx, and more than 90% are of squamous cell 
histologic type. HNC is most often associated with extensive 
lifetime exposure to tobacco and alcohol consumption, although 



high-risk human papillomavirus (HPV) has recently emerged as 
the primary etiologic factor responsible for an increasing number 
of oropharyngeal tumors. ^'^ 

Epidemiologic evidence also supports the hypothesis that 
diet modulates risk, progression and prognosis of head and neck 
cancer.^'^ However, the molecular mechanisms by which dietary 
compounds exert their effects are not entirely understood. 
Epigenetic dysregulation is a key mechanism in tumorigenesis 
that may be influenced by dietary intake. Aberrant DNA meth- 
ylation in cancer is typically characterized by gene promoter 
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Table 1. Clinical characteristics of the study participants (n = 49) 



Patient 
characteristic 




N (%) 


Moan l^r\\ 

iviean \ou)f 
median (range) 


Age 






(41-75) 


Gender 


Male 


39 (80%) 




Female 


10(20%) 






1 


1 (2%) 




Stage 


2 


6(12%) 




3 


12(24%) 






4 


30 (61%) 






Oral cavity 


9(18%) 




Cancer site of first 
primary 


Oropharynx 


26 (53%) 




Hypopharynx 


2 (4%) 




Larynx 


10(20%) 






Other 


2 (4%) 




Tumor tissue HPV 
(+) status 




21 (43%) 






Never 


9(18%) 




Smoking status 


Past 


32 (65%) 






Current 


8(16%) 






AUDIT > 8 and 






Alcohol problem 


drank within 
1 year 


16(32%) 





Specific hypermethylation, leading to gene silencing and hypo- 
methylation of intergenic regions, leading to overall genomic 
instability.^'^ 

Dietary intake has been widely hypothesized to influence 
epigenetic state by determining the availability of functional 
groups involved in the covalent modification of DNA and his- 
tone proteins. ^^'^^ Specifically, nutrients that act as methyl donors, 
such as methionine and folate, or nutrients that act as cofactors 
in one carbon metabolism pathway, such as betaine, vitamin 
B12 and choline, are required for establishment and mainte- 
nance of methylation marks both on DNA and histone protein 
tails. ^^'^^ A deficiency or excess of any of the nutrients involved 
in one-carbon metabolism can potentially alter the availability 
of S-adenosylmethione (SAM), the universal donor involved in 
DNA methylation, and ultimately alter epigenetic state. ^ In addi- 
tion to diet, it has been shown that in smokers, chronic DNA 
damage coupled with reduced DNA repair capacity may induce 
gene promoter methylation.^^ Given these findings, we hypoth- 
esize that, in addition to micronutrients involved in one-carbon 
metabolism, micronutrients involved in cell growth and differ- 
entiation, such as vitamin A, or with antioxidant activity such as 
vitamin C, vitamin E and |3-carotene may indirectly influence 
DNA methylation by reducing DNA damage. 

The objective of this study was to test the hypothesis that pre- 
treatment dietary intake of methyl donors, antioxidants and foods 
abundant in these micronutrients is associated with tumor DNA 
methylation in head and neck cancer. A better understanding of 



how dietary intake is associated with tumor DNA methylation 
in head and neck cancer will enhance our knowledge of disease 
pathogenesis and may allow for the development of individual- 
ized medical nutrition therapy regimens. 

Results 

Descriptive statistics: study sample. The mean age of the 49 
study participants with both dietary intake and methylation data 
was 57.3 y (range 41-75 y); 80% of the participants were male. 
The majority of the cancer sites were oropharynx (53%), evi- 
denced by the relatively high proportion of HPV(+) cases (43%). 
Approximately 85% of the study participants reported being cur- 
rent or former smokers, with an estimated one-third reporting an 
alcohol problem. Clinical characteristics are displayed in Table 1. 

Dietary intake associations with tumor suppressor methyla- 
tion. Reported average daily intakes of micronutrients and food 
group servings reported by study participants are presented in 
Table 2. The distribution of the tumor suppressor methylation 
score is shown in Figure SI. Table 3 shows associations across 
quartiles of micronutrient and food group intake and tumor sup- 
pressor DNA methylation. Higher antioxidant, methyl donor 
and food group intake tended to be negatively associated with 
tumor suppressor score, the sum of individual hypermethylated 
tumor suppressor associated CpG sites, suggesting an inverse 
relationship between DNA methylation at promoter regions of 
these genes and micronutrient intake. The regression coefficients 
in Table 3 can be interpreted as the number of CpG sites of genes 
associated with tumor suppression activity that are methylated 
more than 50% compared with individuals in the lowest quar- 
tile of intake. Specifically, individuals in the lowest quartiles of 
dietary folate, vitamin B12 and vitamin A intake were found to 
have significantly higher tumor suppressor methylation scores 
compared with individuals in the highest quartiles of intake. 
Additionally, tests for trend across quartiles of intake for dietary 
folate, vitamin B12 and vitamin A intake indicated statistically 
significantly higher methylation score with each decreasing 
quartile (Table 3). Individuals in the lowest quartile of crucifer- 
ous vegetable intake had a significantly higher tumor suppres- 
sor methylation score than individuals in the highest quartile of 
intake, and the test for trend across quartiles was also statistically 
significant (p = 0.013). Interestingly, individuals consuming the 
highest servings per day of refined grains had significantly higher 
tumor suppressor methylation score than individuals consuming 
the lowest servings per day of refined grains. However the test for 
trend across quartiles of intake for refined grains was not statisti- 
cally significant. 

In a sensitivity analysis, in which we categorized micronutri- 
ent intake into tertiles as opposed to quartiles and adjusted for 
the same covariates, our findings were similar to the findings 
reported above in that individuals in the lowest tertile of vita- 
min B12, vitamin A and cruciferous vegetable intake had sig- 
nificantly higher tumor suppressor methylation score than those 
in the highest tertile of intake (p = 18.4, (3 = 18.2 and (3 = 18.2, 
respectively; data not shown). The tests for trend across tertiles 
of intake for the relationship between vitamin B12, vitamin A 
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Table 2. Summary of mean daily micronutrient and food group servings by quartile of consumption 



1 


Dietary component 


Quartile 1 


Quartile 2 


Quartile 3 


Quartile 4 






n = 12 


n = 12 


n = 13 


n = 12 


Micronutrients (energy-adjusted) 




Total folate (iJig/d) 


325 


441 


706 


1294 




Food folate (iJig/d) 


202 


269 


330 


422 




Vitamin A (lU/d) 


3869 


6705 


10594 


222111 




Vitamin B6 (mg/d) 


1.8 


2.5 


3.8 


10.5 




Vitamin B12(|jLg/d) 


4.7 


8.5 


13.1 


32.2 




Vitamin C mg/day 


64 


114 


243 


1086 




Vitamin D(IU/d) 


77 


183 


462 


835 




Vitamin E (mg/d) 


5.4 


8.9 


24.6 


288 




Beta-carotene {[xg/d] 


1022 


1913 


3353 


9178 




Methionine (g/d) 


1.4 


1.8 


2.3 


2.8 


Food groups (servings/d) 




Cruciferous vegetables 


0.02 


0.1 


0.22 


0.62 




Green leafy vegetables 


0.07 


0.15 


0.3 


0.94 




Dark yellow & orange vegetables 


0.03 


0.09 


0.19 


0.45 




Refined grains 


0.32 


0.67 


1.17 


2.9 




Whole grains 


0.07 


0.36 


0.86 


2.11 




Red meat 


0.18 


0.4 


0.66 


1.43 


1 


Legumes 


0.07 


0.27 


0.27 


0.71 



and cruciferous vegetables and tumor suppressor methylation 
score were also significant (P , = 0.015, P , = 0.034 and P , = 

C ^ trend trend trend 

0.022, respectively). The inverse association between tumor sup- 
pressor methylation score and food folate was no longer statisti- 
cally significant (p = -14.7 for tertile 3 vs. tertile 1, P^^^^^ = 0.074). 

CpG site-specific associations with vitamin A and vitamin 
B12 intake. With respect to vitamin A intake, there were no 
significant associations with site-specific DNA methylation for 
HPV(+) individuals (all adjusted p values > 0.15). In contrast, a 
number of CpG sites were found to be differentially methylated 
for HPV(-) individuals with the highest levels of both vitamin 
A and vitamin B12 intake compared with those with the low- 
est levels of intake. Interestingly, vitamin B12 intake was also 
found to be associated with hypomethylation of a CpG site asso- 
ciated with RARE (retinoic acid receptor (3) in HPV(+) individu- 
als after adjusting for multiple comparisons (adjusted p value 
< 0.05). The ten sites identified as the most significantly differ- 
entially methylated by both vitamin A and vitamin B12 intake in 
HPV(-) individuals, as well as the one site differentially methyl- 
ated by vitamin B12 intake in HPV(+) individuals are shown in 
Table 4. Consistent with the findings of the methylation score 
analysis, the most differentially methylated CpG sites were, by 
and large, found to be hypomethylated in tumors from indi- 
viduals in the highest quartile of vitamin A and vitamin B12 
intake compared with tumors from individuals in the lowest 
quartiles of intake. There was one exception, CHGA (chro- 
mogranin A), which was found to be hypermethylated in 
tumors from HPV(-) negative individuals who consumed the 
largest amounts of vitamin B12. CpG sites associated with 



C4B (complement component 4B), GML (glycosylphosphati- 
dylinositol anchored molecule like protein), MLFl (myeloid leu- 
kemia factor 1), and THPO (thrombopoietin) were found to be 
significantly associated with methylation in tumors from HPV(-) 
individuals based on both vitamin A and vitamin B12 intake. 

Table SI presents the results of the ranked gene set enrich- 
ment analysis (GSEA). An analysis of positional enrichment of 
genes differentially methylated by vitamin B12 intake in HPV(+) 
cancers revealed a trend toward hypomethylation of chromosome 
12pl2 with increased vitamin B12 intake. The same chromo- 
somal region was identified as hypomethylated in tumors from 
HPV(-) individuals with the highest levels of vitamin A intake, as 
was chromosome 8p21. While there was no consistent enrichment 
of gene ontology biological processes or Kyoto Encyclopedia of 
Genes and Genomes (KEGG) pathways based on dietary intake 
for HPV(-) tumors, for HPV(+) tumors there was a significant 
enrichment of processes involved in immune function based on 
vitamin B12 intake. 

Discussion 

To our knowledge, this is the first study to comprehensively 
examine the association between dietary intake and promoter 
methylation of genes related to head and neck cancer. These 
novel findings suggest diet is significantly associated with epi- 
genetic events that occur in head and neck cancer. While previ- 
ous studies have indicated deficiencies of folate, vitamin B12 and 
vitamin A are associated with increased risk of developing head 
and neck cancer^'^^'^^ as well as with poorer prognosis in patients 
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Table 3. Results of unadjusted and adjusted linear regression showing associations between tumor suppressor methylation score and quartile of 
energy adjusted micronutrient intake or intake of food group, setting quartile 1 (lowest intake) as the reference category 



Unadjusted 




Micronutrient 


Quartile 2^ (S.E) 


Quartile 3^ (S.E) 


Quartile 4^ (S.E) 


p value for trend 




Dietary folate 


-4.3 (8.6) 


-7.0 (8.6) 


-16.8 (8.4) 


0.047 




Total folate 


1.3 (8.7) 


-1.0 (8.7) 


-6.2 (8.7) 


0.450 




Vitamin A 


-8.6 (8.6) 


-7.9 (8.6) 


-13.3 (8.6) 


0.146 




Vitamin B6 


-7.1 (8.6) 


-1.0 (8.6) 


-11.4 (8.6) 


0.301 




Vitamin B12 


-7.2 (8.0) 


-6.6 (8.2) 


-23.9 (8.2) 


0.008 




Vitamin C 


0.8 (8.6) 


2.0 (8.6) 


-8.4 (8.6) 


0.389 




Vitamin D 


-8.8 (8.6) 


-11.2 (8.6) 


-10.2 (8.6) 


0.218 




Vitamin E 


-1.8(8.8) 


-3.9 (8.8) 


-2.4 (8.8) 


0.727 




Beta-carotene 


-3.3 (8.6) 


-7.5 (8.6) 


-13.1 (8.6) 


0.111 




Methionine 


1.5 (8.5) 


11.7 (8.5) 


-4.4 (8.5) 


0.932 




Cruciferous vegetables 


-4.7 (8.0) 


-16.5 (8.0) 


-18.1 (8.5) 


0.014 




Green leafy vegetables 


0.7 (8.5) 


-6.1 (8.9) 


-6.6 (8.9) 


0.343 




Dark yellow vegetables 


-5.3 (9.3) 


-5. 3 (9.0) 


-10.1 (8.4) 


0.245 




Refined grains 


12.3 (8.5) 


9.3 (8.8) 


14.3 (8.7) 


0.150 




Whole grains 


9.1 (8.7) 


8.0 (9.3) 


-1.5 (8.8) 


0.808 




Red meat 


-17.8 (8.2) 


-11.3 (8.5) 


-12.4 (8.0) 


0.185 




Legumes 


-5.6 (8.3) 


-5.4 (9.6) 


-9.8 (9.0) 


0.307 


Adjusted^ 




Dietary folate 


-3.1 (8.9) 


-4.0 (9.5) 


-18.7 (8.9) 


0.043 




Total folate 


4.6 (9.1) 


-3.2 (9.6) 


-9.1 (9.2) 


0.258 




Vitamin A 


-13.3 (9.8) 


-15.1 (9.9) 


-21.1 (9.7) 


0.040 




Vitamin B6 


-5.5 (9.6) 


-2.2 (9.3) 


-12.9 (9.2) 


0.219 




Vitamin B12 


-7.6 (8.8) 


-11.0 (8.7) 


-27.6 (8.8) 


0.003 




Vitamin C 


0.2 (8.9) 


-1.9 (9.4) 


-16.6 (9.3) 


0.089 




Vitamin D 


-11.4 (9.3) 


-15.6 (9.5) 


-18.1 (9.8) 


0.059 




Vitamin E 


-4.3 (9.2) 


-11.4 (9.8) 


-12.2 (9.9) 


0.175 




Beta-carotene 


-3.4 (9.2) 


-12.5 (9.0) 


-15.6 (9.4) 


0.055 




Methionine 


3.9 (9.1) 


9.6 (9.1) 


-7.9 (9.7) 


0.687 


r 


Cruciferous vegetables 


-6.8 (8.5) 


-16.2 (8.5) 


-20.4 (9.2) 


0.013 




Green Leafy vegetables 


4.8 (9.0) 


-8.7 (9.4) 


-10.9 (9.5) 


0.183 




Dark Yellow vegetables 


-6.8 (9.6) 


-15.4 (10.8) 


-15.4 (9.8) 


0.098 




Refined grains 


15.3 (9.7) 


13.5 (10.3) 


22.6(10.7) 


0.069 


[ 


Whole grains 


10.2 (9.7) 


0.8 (10.7) 


-3.8 (10.9) 


0.542 




Red meat 


-14.8 (10.3) 


-12.8 (10.5) 


-14.9 (11.8) 


0.189 




Legumes 


-6.2 (8.9) 


-4.5 (10.6) 


-14.2(10.6) 


0.227 



^Adjusted for age, sex, smoking status, HPV status and problem drinking. Food group analyses also adjusted for total calorie consumption. Bold text 
indicates p < 0.05. 



diagnosed with the disease, our study has shown a decrease in 
tumor suppressor methylation in the highest quartile of vitamin 
B12 and vitamin A intake. Epigenetic dysregulation, such as the 
silencing of tumor suppressor genes via promoter hypermethyl- 
ation, is a major event in malignant disease. 

Though the underlying biological mechanism related to how 
folate and vitamin B12 deficiencies may be associated with DNA 
hypermethylation has not been elucidated, we have two suggested 



hypotheses. First, folate and vitamin B12 play important roles in 
the synthesis of DNA and DNA damage repair.^^'^^ Increases in 
DNA methylation have been observed in an experimental model 
inducing in vitro DNA damage, in the form of DNA strand 
breaks, and homology directed DNA repair.^^ An increased rate 
of DNA strand breaks was also observed in both animal models 
and cell culture models of folate deficiency,^^'^^ potentially pro- 
viding a mechanism that links folate deficiency, DNA strand 
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Table 4. CpG sites with tumor DNA metiiyiation values significantly different between highest (4th quartile) and lowest (1st quartile) of vitamin B12 and 
vitamin A intake, stratified by HPV status 



Gene symbol 


Chromosome 


CpG coordinate 


Distance to TSS t-value 


p value 


Adjusted p value 


HPV positive cases: methylation by vitamin B12 intake 


RARB 


3 


25444872 


114 


-5.72 


3.18E-05 


0.0479 


HPV negative Ceases: methylation by vitamin B12 intake 


1 C4B 


6 


32057984 


171 


-5.39 


1.79E-05 


0.0178 


GML 


8 


143912938 


-281 


-5.23 


2.64E-05 


0.0178 


PTHR1 


3 


46893982 


-258 


-5.11 


3.55E-05 


0.0178 


MLFl 


3 


159771921 


243 


-4.78 


8.14E-05 


0.0306 


THPO 


3 


185578143 


483 


-4.58 


0.0001 


0.0341 


CHGA 


14 


92459297 


52 


4.57 


0.0001 


0.0341 


XIST 


X 


72989344 


-31 


-4.44 


0.0002 


0.0350 


LCN2 


9 


129951398 


-141 


-4.41 


0.0002 


0.0350 


SRC 


20 


35406338 


-164 


-4.40 


0.0002 


0.0350 


SFN 


1 


27062338 


118 


-4.26 


0.0003 


0.0442 


HPV negative cases: methylation by vitamin A intake 


PTHLH 


12 


28016940 


-757 


-6.55 


1.11E-06 


0.0017 


LM02 


11 


33870264 


148 


-5.80 


6.64E-06 


0.0050 


SIN3B 


19 


16800704 


-514 


-5.26 


2.48E-05 


0.0092 


TNFSF8 


9 


116732775 


-184 


-5.19 


2.91 E-05 


0.0092 


C4B 


6 


32057984 


171 


-5.17 


3.06E-05 


0.0092 


GML 


8 


143912938 


-281 


-4.66 


0.0001 


0.0241 


HGF 


7 


81238681 


-1293 


-4.65 


0.0001 


0.0241 


MLFl 


3 


159771921 


243 


-4.59 


0.0001 


0.0244 


THPO 


3 


185578143 


483 


-4.51 


0.0002 


0.0244 


HOXAn 


7 


27191320 


35 


-4.44 


0.0002 


0.0244 



breaks, and DNA methylation. Second, while the DNA methyla- 
tion changes that occur in cancer have been well characterized, 
there is also interest in describing the DNA methylation changes 
that occur with aging, as these changes could reflect what occurs 
in the early stages of carcinogenesis. A recent study has described 
that DNA methylation of two genes, ERa and SFRPl, in colon 
tissue increases with age and also with folic acid supplementation 
status. In our study, we not only observed a decrease in tumor 
suppressor methylation in the highest quartile of vitamin A, vita- 
min B12 and food folate intake, but also a trend toward increased 
tumor suppressor methylation in the individuals who consumed 
the highest levels of refined grains. During the refinement process 
these grains are fortified with folic acid, niacin, riboflavin, thia- 
min and iron. It is unclear how these nutrients may be associated 
with increased DNA methylation and more research is necessary. 
The relative bioavailability of these nutrients could modulate the 
normal stochastic changes in methylation that occur due to aging 
and may be reflected in tumor DNA methylation. 

Other studies support our finding that increased consump- 
tion of folate and folate-rich foods are associated with decreased 
methylation of genes associated with tumor suppression. Another 
study conducted in head and neck cancer patients reported 
increased methylation of pW^^^"^, a gene involved in cell cycle 



regulation that is commonly lost or silenced in head and neck 
cancer via epigenetic mechanisms, in those consuming the lowest 
levels of folate compared with the highest levels. 

Based on our findings that folate intake was inversely associ- 
ated with tumor suppressor DNA methylation, it is not surprising 
that high cruciferous vegetable consumption was associated with 
lower methylation levels, given this food group is a rich source of 
folate. However, it is interesting to note that cruciferous vegeta- 
bles are also rich sources of glucosinolates such as sulforphane and 
indole-3-carbinol, both of which are believed to exhibit anti-can- 
cer activities. Though there is little evidence suggesting involve- 
ment of sulforaphane in DNA methylation, this compound has 
been shown to act as a histone deacetylase (HDAC) inhibitor. 
HDACs have been shown to be widely overexpressed in a number 
of cancers and are involved in gene silencing by inducing a more 
condensed chromatin state, leaving DNA less accessible for tran- 
scription.^^ HDAC inhibitors, such as sulforaphane, have been 
shown to be promising cancer therapies because of their ability to 
induce expression of epigenetically silenced genes. 

The finding that vitamin A intake is significantly associ- 
ated with decreased DNA methylation in HPV(-) tumors is 
unexpected, considering vitamin A does not play a direct role 
in the one-carbon metabolism pathway. However, it is possible 
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vitamin A can alter DNA methylation by influencing the avail- 
ability and activity of methyltransferases. Two studies have 
reported rats treated with all-trans-retinoic acid (ATRA), a 
vitamin A metabolite, displayed increased levels and activ- 
ity of glycine A^-methyltransferase (GNMT), the enzyme 
involved in the conversion of S-adenosyl methionine (SAM) to 
S-adenosylhomocysteine (SAH).^^'^^ Increased GNMT activity 
led to a reduction of SAM, thus potentially resulting in a reduc- 
tion of methylation. Similarly, Das et al. reported demethylation 
of 402 gene promoters in neuroblastoma cells treated with ATRA. 
These cells also exhibited downregulation of the methyltransfer- 
ases, DNMTl and DNMT3B, which the authors proposed as the 
potential mechanism by which ATRA exerts its effect on DNA 
methylation.^^ 

The GSEA findings of this study confirm the previous work 
of our group and others describing the differences in tumor biol- 
ogy between HPV(+) and HPV(-) HNCs.^^'^^ In this study, we 
identified that gene sets involved in immune function, including 
cytokine production and antigen processing and presentation, 
were significantly enriched based on DNA methylation differ- 
ences in HPV(+) tumors from individuals with high and low 
levels of vitamin B12 intake. There was no consistent KEGG or 
Gene Ontology Biological Processes enriched in HPV(-) tumors, 
based on either vitamin A or vitamin B12 intake. A positional 
enrichment analysis identified the same chromosomal region, 
chromosome 12pl2, as differentially methylated in HPV(+) 
tumors based on vitamin B12 intake as HPV(-) tumors based on 
vitamin A intake. This region has previously been identified as 
a candidate tumor suppressor region in head and neck and other 
cancers, and may represent a region of the genome that is more 
labile to the epigenetic influences of dietary intake. 

Results of this preliminary study should be interpreted in 
light of several limitations. Our small sample size poses analytic 
limitations in detecting associations of interest and the complex- 
ity of statistical tests used. Although we have adjusted for several 
covariates, there are no doubt other unmeasured sources of con- 
founding we were not able to consider. The FFQis susceptible to 
measurement error and potential unquantified systematic biases 
that may have led to misclassification of dietary intake. The FFQ 
was designed to assess usual dietary intake over the course of 
the year prior to diagnosis with HNC. Thus we cannot make 
any conclusions about the association between lifetime dietary 
exposures and tumor DNA methylation. While the Goldengate 
Cancer Panel provides a comprehensive measurement of meth- 
ylation across 1,505 CpG sites in promoter regions of genes with 
known roles in cancer, this technology does not provide an unbi- 
ased assessment of the epigenome. Finally, due to these limita- 
tions, the findings of this study should be validated in another 
set of tumors from similar individuals. 

In conclusion, dietary intake of micronutrients involved in 
one-carbon metabolism such as folate and vitamin B12, as well 
as other micronutrients such as vitamin A, may be significant 
regulators of epigenetic events occurring in head and neck can- 
cer. These findings have potential clinical implications for the 
treatment of head and neck cancer. The idea that dietary inter- 
ventions could potentially reprogram the epigenome in such a 



way as to optimize the likelihood of positive disease outcomes 
is appealing. Further research is necessary to better understand 
the complex role these micronutrients play in regulating the epi- 
genetic process, as well as how nutrients may differentially regu- 
late this process according to HPV status. Such studies could 
give rise to highly specific randomized controlled trials and ulti- 
mately to the development of individualized medical nutrition 
therapy regimens that may improve prognosis in the head and 
neck cancer population. 

Methods 

Design. This was a cross-sectional study of patients enrolled in 
the University of Michigan Head and Neck Cancer Specialized 
Program of Research Excellence (SPORE). The independent 
variables were micronutrients involved in one-carbon metabo- 
lism, antioxidant nutrients and food groups abundant in these 
micronutrients. Confounding variables included age, sex, smok- 
ing, alcohol problem and HPV status. The outcome variable 
was DNA methylation. Our analyses were limited to foods and 
nutrients we hypothesized are associated with DNA methyla- 
tion — namely those involved in one-carbon metabolism or with 
antioxidant activity. A two-tiered analysis was conducted; first 
we observed whether dietary intake is associated with methyla- 
tion at CpG sites in tumor suppressor genes overall and second, 
we observed whether DNA methylation at specific CpG sites is 
associated with dietary intake. 

Study population. Sixty-eight newly diagnosed cases of head 
and neck cancer were identified and recruited at two sites partici- 
pating in the University of Michigan Head and Neck SPORE: 
the University of Michigan Health System and the Ann Arbor 
Veterans Administration (VA) Medical Center. Individuals eli- 
gible for participation included patients diagnosed with first 
primary head and neck cancer between January 1, 2003 and 
December 31, 2005 who signed written, informed consent, 
completed an epidemiologic questionnaire that includes ques- 
tions on demographic and behavioral data, and had a paraffin- 
embedded tumor available for analysis with adequate residual 
tissue for microdissection. Exclusion criteria included: (1) < 18 
y of age, (2) pregnant, (3) non-English speaking, (4) diagnosed 
as mentally unstable or (5) a diagnosis of another non-upper 
aerodigestive tract cancers (such as thyroid or skin cancer). For 
the purposes of the current analysis, patients were excluded if 
they had not completed a baseline FFQ (n = 19), yielding a final 
sample size of 49 participants. Excluded patients did not differ 
significantly from included patients in respect to age, sex, tumor 
site, cancer stage or smoking status. This study was approved as 
being within the ethical standards of the Institutional Review 
Boards of the University of Michigan and the Ann Arbor VA 
Medical Center. 

Tumor blocks were recut for uniform histopathologic review 
and microdissection, with the first and last slides of a series of 
12 reviewed by a qualified pathologist (JM and ND) to confirm 
the original diagnosis and to circle areas for DNA extraction. 
Percent cellularity was estimated for each tumor, and areas with 
> 70% cancer cellularity were designated for use in the analyses. 
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Measures. FFPE tissue, DNA isolation and bisulfite conver- 
sion. Regions that were identified for DNA extraction were cored 
from the formalin fixed paraffin embedded (FFPE) tissue blocks 
using an 18 gauge needle. Isolation of DNA from the cored tis- 
sue samples was performed using the QIAamp DNA FFPE 
Tissue Kit (Qiagen) with a modification to the manufacturer's 
recommended lysis protocol (incubation overnight at 56°C in 
lysis buffer) . DNA concentration and purity was confirmed via 
NanoDrop spectrophotometer (Thermo Scientific). We per- 
formed sodium bisulfite modification on 500 ng to 1 |JLg of DNA 
using the EZ DNA Methylation kit (Zymo Research) following 
the manufacturer's recommended protocol. 

Bead array methods. The commercially available lUumina 
Goldengate® Methylation Cancer Panel was utilized to detect 
DNA methylation patterns in tumor samples. The Cancer Panel 
interrogates 1,505 CpG sites located in known CpG islands 
across 807 genes related to cancer, including oncogenes, tumor 
suppressor genes, imprinted genes, and genes involved in cell 
cycle regulation, DNA repair, apoptosis and metastasis. Bead 
arrays were run at the University of Michigan DNA Sequencing 
Core Facility according to the manufacturer's protocol. Briefly, 
bisulfite converted tumor DNA was hybridized to the bead array 
as described previously,^^ and bead arrays were imaged using 
Illumina Cancer Panel Reader software. Raw bead array fluo- 
rescence data was initially analyzed using Illumina BeadStudio 
Methylation software, which converts fluorescence values of the 
methylated (Cy5) and unmethylated (Cy3) alleles into an aver- 
age methylation value at a specific probe using the formula (3 = 
[Max(Cy5,0)]/[Max(Cy5,0) + Max(Cy3,0) + 100] with pe(0,l). 

Methylation at specific CpG probes on the Goldengate Cancer 
Panel has been shown to be biased by probe thermodynamic 
properties.^^ Known biases include probe length and GC con- 
tent, which can affect the melting temperature of the probes as 
well as probe fluorescence intensities. Thus, we used the method 
proposed by Kuan et al. to normalize our average (3 values based 
on probe length and GC content.^^ 

Detection p values on the Goldengate Cancer Panel are cal- 
culated based on fluorescence signal at a probe compared with 
background fluorescence and represent one minus the prob- 
ability) that a signal is stronger than background fluorescence. 
The weighted methodology proposed by Kuan et al. was used to 
develop sample and site weights based on p-values of detection. 
Both samples and sites with larger detection p-values are gener- 
ally considered less reliable and were down-weighted in further 
gene specific analyses. 

Dietary intake estimation. Dietary components of interest in 
this study included antioxidant micronutrients, micronutrients 
involved in one-carbon metabolism, and food groups serving as 
rich sources of these micronutrients including cruciferous vegeta- 
bles, green leafy vegetables, dark yellow vegetables, refined grains, 
whole grains, red meat and legumes. Pre-treatment dietary intake 
was collected using a self- administered, semi-quantitative FFQ, 
designed to assess respondents' usual dietary intake from food and 
supplements over the year prior to diagnosis. The reproducibility 
and validity of this 131-item questionnaire has been reported 
previously.^^'^^ The FFQ includes given standard portion sizes 



for each item (i.e., 1 orange or 3-5 oz chicken) which allowed 
participants to choose their average frequency of consumption 
over the past year from a list of 1-9 choices ranging from "almost 
never" to "> 6 times per day." Daily energy and nutrient intake 
was estimated by summing intakes from each food based on the 
given standard portion size, reported frequency of consumption, 
and nutrient content of each food item.^^ Daily food group serv- 
ings were estimated by summing the frequency weights of each 
food item based on reported daily frequencies of consumption. 

Covariates. Each participant completed a self-administered 
health questionnaire at baseline. This questionnaire was designed 
to collect demographic data as well as data on smoking, alcohol 
use, physical activity, sleep, comorbidities, depression and qual- 
ity of life. Age, sex, smoking and alcohol data were collected 
from the health questionnaire. Smoking data permitted analy- 
sis by "current-past-never" smoking. Alcohol consumption was 
measured using the previously validated Alcohol Use Disorders 
Identification Test (AUDIT) .^^ An AUDIT score > 8 was con- 
sidered problem drinking. Tumor site and stage was recorded 
from operative notes and surgical pathology forms at baseline. 
As we have previously shown HPV status to be associated with 
both tumor DNA methylation profiles at specific gene promoter 
regions^^ and micronutrient intake,^*^ we considered HPV status 
as a potential confounder in statistical analyses. 

HPV status was determined by an ultra-sensitive method 
using real-time competitive polymerase chain reaction and 
matrix-assisted laser desorption/ionization-time of flight mass 
spectroscopy with separation of products on a matrix-loaded sili- 
con chip array. Multiplex PGR amplification of the E6 region of 
15 discrete high-risk HPV types (HPV 16, 18, 31, 33, 35, 39, 45, 
51, 52, 56, 58, 59, 66, 68 and 73), and human GAPDH control 
was run to saturation followed by shrimp alkaline phosphatase 
quenching. Amplification reactions included a competitor oligo 
identical to each natural amplicon except for a single nucleo- 
tide difference. Probes that identify unique sequences in the 
oncogenic E6 region of each type were used in multiplex single 
base extension reactions extending at the single base difference 
between wild-type and competitor HPV so that each HPV type 
and its competitor were distinguished by mass when analyzed on 
the MALDI-TOF mass spectrometer. 

Statistical methods. Methylation score. To create a measure- 
ment of stochastic methylation across genes that could be asso- 
ciated with dietary intake, we developed a novel method of 
calculating a methylation score based on methylation of genes 
associated with specific cellular functions. Since the Goldengate 
Cancer Panel quantifies methylation of cancer-related genes with 
various functions, we first subsetted genes associated with our 
main cellular function of interest, tumor suppression. To identify 
genes associated with tumor suppression, we conducted a sim- 
ply query of the National Center for Biotechnology Information 
(NCBI) Gene database for "tumor suppressor" and extracted 
Gene IDs. A total of 444 individual CpG sites across 237 genes 
on the Goldengate Cancer Panel were linked to tumor suppres- 
sion. Next, for each study participant, the methylation score was 
defined as the number of tumor suppressor gene-associated CpG 
sites with normalized average (3 values above 0.5. A (3 value cutoff 
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of 0.5 was established as a methylation level likely to result in 
decreased expression of tumor suppressor genes. Thus, the meth- 
ylation score ranged from 0 to 444. We infer that a higher meth- 
ylation score is associated with higher levels of DNA methylation 
of promoter regions of genes associated with tumor suppression 
and cell cycle regulation. 

Micronutrients were adjusted for total energy intake using 
the residual method.^*^ Descriptive statistics, including means, 
medians, maxima and minima were calculated for both the 
energy adjusted nutrients and the methylation score. Associations 
between energy- adjusted micronutrient intake, food group intake, 
and methylation score were assessed with Spearman rank correla- 
tion due to the typical non-normal distribution of micronutrient 
and food group intake. Additionally, we used linear regression 
models with methylation score as the outcome and energy- 
adjusted micronutrient or food group intake as the main predic- 
tor, categorizing dietary variables into quartiles, and adjusting 
for age, sex, smoking status, problem alcohol drinking, and HPV 
status, as well as analyzing for trend. To examine the robustness 
of our findings, we performed a sensitivity analysis in which we 
categorized micronutrient intake into tertiles as opposed to quar- 
tiles and adjusted for the same covariates. Regression diagnostic 
plots were examined to ensure that associations observed were not 
being driven by influential points. All analyses were performed in 
SAS software. Version 9.2 or R Version 2.13.0. 

Site-specific analyses. To determine whether methylation at 
specific sites of the genome may be labile to effects due to dietary 
intake, we calculated associations between DNA methylation at 
specific CpG sites and dietary intake, stratifying by HPV sta- 
tus. To limit the number of multiple comparisons conducted, we 
restricted our analysis to micronutrients identified as significant 
in the tumor suppressor gene methylation score analysis follow- 
ing sensitivity analysis, vitamin B12 and vitamin A, and com- 
pared site specific DNA methylation between highest and lowest 
quartile of intake. 

The association between micronutrient intakes and individual 
CpG site DNA methylation at the 1505 CpG sites measured on 
the Goldengate Cancer Panel were examined using the Limma 
package in R 2.10.1.^^ Sample weights based on detection p-values 
across samples were used in the ImFit function from the Limma 
package to downweight samples with higher detection p-values. 
An empirical Bayes method (using the eBayes function in Limma) 
was used to shrink standard errors to a common value and to 
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